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INTRODUCTION 



One of the most fascinating manifestations of many- 
body physics in mesoscopic systems is the Kondo effect. 
It was originally observed in the 1930s in gold wires and 
formally explained in 1964 by Jun KondcWin metals that 
contain magnetic impurities. The Kondo effect became 
an intense topic of research in the context of mesoscopic 
systems^ after its observation in quantum dots (QDs) in 
the late nineties P The Kondo effect in QDs is revealed 
as a zero-bias conductance enhancement at low temper- 
atures when the QD is spin polarized and strongly cou- 
pled to the conduction electrons (such as a QD that is 
tunnel-coupled to leads). Since its discovery in QDs, the 
Kondo effect has been observed and studied in a variety 
of exciting mesoscopic systems: individual atoms car- 
bon nanotubes, 5 molecules}^ bucky balls^and grapheneP 

Since the Kondo transport can strongly influence the 
Fano phenomenon, there has recently been a lot of in- 
terest in the interplay between these two effects in QDs. 
This interest comes from the use of Fano resonances as 
a tool to study coherence and to provide valuable trans- 
port information}^^ that may result in future nano-scale 
device applications.^^ The Fano effect results from the 
interference between a scattering event within a contin- 
uum of states (background process) and another one by 
means of a discrete energy level (resonant process). A 
Fano resonance manifests itself as an asymmetric conduc- 
tance line-shape!^ Fano resonances have been observed 
in many different physical systems and extensively stud- 
ied in nanoscale systems!^ Scientific awareness of the 
interaction between the Fano and the Kondo effect be- 
gan with STM measurements on metal surfaces,^ how- 
ever, it found its greatest impact in mesoscopic systems. 
The Fano-Kondo effect is observed when Kondo trans- 
port takes over the background scattering process in the 
Fano effect. The trademark of such interaction is a con- 
ductance valley when tuning the energy level of the res- 
onant process. 



So far, only a small number of dedicated Fano-Kondo 
studies have been reportedP^HH] Sasaki et al. 20 per- 
formed transport measurements through a side-coupled 
double QD configuration and found a strong modulation 
of the Kondo conductance by the Fano effect. Their re- 
sults are in good agreement with a tight-binding model 
presented in a parallel theoretical publication.^ In a the- 
oretical work, Fang and Luo^ have calculated the con- 
ductance in an Aharonov-Bohm interferometer contain- 
ing a Kondo QD coupled to a non-interacting QD. They 
predict a Fano-Kondo effect that could decoupled the two 
dots through the manipulation of the Kondo resonance 
by varying the inter-dot coupling. 

In this publication we investigate the Fano-Kondo in- 
terplay in an Aharonov-Bohm ring, containing a non- 
interacting QD and a Coulomb interacting QD. Using 
a scattering matrix approach combined with the slave- 
boson mean-field approximation (SBMF), we calculate 
the conductance and the Fano parameter. The deriva- 
tion of the Fano form for the conductance allows us to 
analyze in-depth the Fano-Kondo phenomenon and, spe- 
cially, the origin of the transport features. In addition, 
we present a qualitative comparison with the experiment 
performed by Verduijn et alM, and suggest an explana- 
tion of their observations in terms of a Fano-Kondo in- 
terference phenomenon in an Aharonov-Bohm geometry. 

The paper is organized as follows: In section [TT] we in- 
troduce the Anderson Hamiltonian for a non-interacting 
QD and a Coulomb interacting QD connected in paral- 
lel. We diagonalize this Hamiltonian resulting in a set 
of self-consistent equations, which we solve numerically 



to obtain the conductance. In section [HI] we present our 
calculations. These lead us to predict a Fano footprint of 
an underlying Kondo channel and to show that the mag- 
netic flux can open a gap between two resonances and 
alter the electron's path preference for transport. We 



conclude in [IV] with a discussion of applications of our 
results. 
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II. MODEL & METHODOLOGY 

We consider two QDs embedded in an Aharonov-Bohm 
ring which is attached to two ID conducting leads. Fig. 
[I] shows a schematic representation of the nanostructure. 
A magnetic field is used to create a magnetic flux ($) 




FIG. 1. Tight-binding model for two QDs embedded in an 
Aharonov-Bohm ring nanostructure. The QDs are labeled 
with indexes u and d for the upper and lower QD, respectively. 

through the ring, which gives a phase contribution ((/>) 
to the electron's wavefunction. The Hamiltonian of the 
system is given by 

H =eF^2c] a c j , a -t J^(ct a c, +li<r + cx -) 

+ e nf\^fn,cr + Udf^fd^fdlfdl 

Rc\ a f n , a +C.c). (1) 

n,cr 

We assume that each QD has one single energy level la- 
beled e n (with n = u for the upper QD and n = d for 
the QD in the lower path). In the leads each site has a 
single-particle energy at the Fermi energy e F . The hop- 
ping energy between neighboring sites in the leads is to. 
Taking into account the Aharonov-Bohm phase picked up 
by the electron ((/>), the hopping between the ID leads 
and the ring structure is given by t n ^ = t n e ±% i and 

tn,R — tn^ 1 *, where t n ^ and t n ,R (n = u,d), describe 
the hopping of an electron from the upper (u) or lower 
(d) QD to the left (L) or right (R) lead. The operators 
Cj, a (Cj a ) and f n ,a(fn,a) correspond to the annihilation 
(creation) operators in the leads and in the QDs, respec- 
tively. Finally a Coulomb interaction is present in the 
lower QD with strength Ud- 

We calculate the conductance of this nanostructure in 
three regimes, depending on e^, the energy level of the 
lower QD, 1) the Kondo regime (e^ < e Fl Ud > e F — Q 
and < e F — Q, where Yd is the tunnel coupling of 
the lower QD to the leads), 2) the mixed- valence regime 
(e d < e F , U d > e F - e d and T d ~ e F - e d ) and 3) 
the empty orbital regime (e^ > e F ). For simplicity, we 
take Ud — >• oo and use the infinite-U slave boson mean- 
field approximation to diagonalize the Hamiltonian H 
[Eq. |Tj)] .l 23 l 24 l Within this scheme we introduce a bo- 
son creation (annihilation) operator b\(bd) that acts on 
the lower QD. Now the b\ |0) state represents the empty 



state in the lower QD. We then have to create a boson 
when removing an electron from the QD to leave it in the 
empty state, f^ a /d^b^. Since Ud —> oo double oc- 
cupancy is forbidden which is imposed by the constraint 
bdpd + ^2 a fd a fd^ = 1- Within this exact representation, 
the Hamiltonian H of the system [Eq. ([I])] becomes: 

H S b = e F^c ] j a c j)CT - £o^(4,<7 C j+i,^ + c.c.) 

+ E tnfLjn,. + Ub% + E flM* ~ 1) 
n,cr <t 

- ^2(tu,LC-l,afu,(T + t Uy Rc\ a f Uy(T + C.C) 

a 

- Y J { ~ t ^-i,Jd,J d + t d , R clj d} J d + c.c).(2) 

(J 

Here the constraint has been added by means of a La- 
grange multiplier, A^. Since the Hamiltonian Hsb [Eq. 
([2|] remains difficult to solve, we use the mean-field ap- 
proximation were fluctuations around the average value 
of the boson operators in Eq. ([2| are neglected. By sub- 
stituting the boson operators with their real expectation 
values, bd —>< bd >, the slave-boson Hamiltonian [Eq. 
([2|] is replaced by an effective mean-field Hamiltonian, 

+ 2(e d + \ d )tif d + \ d (<b d > 2 -1) 

+ 2e„/^/„ - 2(i„ )L cl 1 /„ + t u , R c\f u + c.c) 

-2<b d > (t d:L clj d + t d , R 4fd + cc). (3) 

The factor of 2 takes into account the contribution from 
both spin directions. We now minimize the ground state 
energy of Eq. ^ with respect to the Lagrange multiplier, 
X d , and the boson expectation value, < b d >, 



d < Hsb,mf > 
d\ d 

d < Hsb,MF > 



d<b d > 



2 < flfd > + < b d > 2 -1 = 0(4a) 
A d < b d > 2 -R d = 0, (4b) 



where 



R d = <b d > (t diL < clj d > +t dtR <c\f d > 

+ t\ L < /]c_! > +tl R < /] Cl >). (5) 

The fermion expectation values in Eq. Q are obtained 

by diagonalizing the mean-field Hamiltonian [Eq. 

We take as ansatz the linear combination of atomic or- 

bitals, \V k ) = Y,n a n \ n ) + Y,j a j \j), where a j and a n 
are the probability amplitudes to find the electron, with 
momentum k and energy w = e F — 2£q cos (/c), at site j in 
the leads and at quantum dot n = u,d, respectively. We 
assume that the transport of electrons is described by an 
incoming plane wave that is reflected and transmitted at 
the Aharonov-Bohm nanostructure. Assuming an incom- 
ing electron from the left lead (L), where the chemical 
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Vl, yields: a k L = e ik ' j 



TLLe 



-ik-j 



potential is at ep 
for j < and a k L = TRL,e lk 'i for j > 0, where < /c < 

= arccos(|^). and r^L correspond to the trans- 
mission amplitude from the left to the right lead and t he 
reflection amplitude in the left lead, respectively! 25 * 26 *. 

By substituting the ansatz, into the Schrodinger 

equation using the mean-field Hamiltonian [Eq. (|3|], we 
arrive at the following system of linear equations: 

- ^o^ll +^e*4a uL + t d e z *a dL =t (6a) 
-toTRL + l 4 a nL + t d e l ±a dL = (6b) 



+(w - e u )e lk a uL 



t d e l *r LL +t d e l ±r RL 



+{w - e d )e-* k a dL = -t d e** e~ M ,(6d) 

where we have renormalized the coupling and the en- 
ergy level of the lower QD as ^ =< b d > t d and 
£d — e d + Ad, respectively. By solving Eq. (p|) we ob- 



tain all the probability amplitudes, a k L and cQ L , for an 
incoming wave from the left. In the same way we solve for 
an incoming wave from the right lead, where the chem- 
ical potential is at ep + Vr, by using the coefficients 
a k R = e~ ik ' 3 ' +r RR e ik ' j for j > and a k R = r LR e~ ik ' j for 

j < 0, where < k < k R = arccos(^). Since the mag- 
netic field breaks the time reversal symmetry, coefficients 
and probability amplitudes for a left and right incoming 
wave are different. With all the coefficients we can com- 
pute the fermion expectation values in Eq. Q. 2 ^ We 
obtain the boson parameter < b d > and Lagrange mul- 
tiplier Ad from the self-consistent set of equations [Eqs. 
(?? -[6])], at equilibrium and zero temperature. From the 
transmission coefficient in Eq. (|6| the conductance at the 
Fermi energy (k = %) reads G = ^- \t R l(%)\ 2 , with 



ground scattering amplitude. When \q\ » 1 the elec- 
tronic transport is dominated by the resonant process 
and the conductance takes a Lorentzian shape. As the 
background process starts to contribute, the Fano param- 
eter increases. When both processes contribute signifi- 
cantly \q\ ~ 1, and the line-shape becomes asymmetric. 
When the background process takes over as the main 
mechanism for transport \q\ ~ and the conductance 
shows a symmetrical dip centered at the resonance.^ 
Below, we show that by tuning the energy level of the 
interacting QD we can set the conductance to any of 
the line-shapes mentioned before, and that interference 
tuned by the magnetic flux leads to interesting behavior 
of the transmission. First we use Eq. ^ to write the 
conductance, in terms of e n , in a Fano formP^ 



G = G Q 



|(6n+e )+gr| 2 

(e n + e ) 2 +T 2 



(8) 



where the amplitude G Q , the resonant point e D , the peak 
broadening T and the Fano parameter q are given by: 



Go 



2e 2 1 
h t? 2 + 1' 



e G = r„ 2 \ - cos (-), 



r = r ,, 2 +sin 2 (f 



T] 2 + 1 

rj[r] 2 cos(</>) — sin 2 (^) — i(r] 2 + 1) sin(0)] 
f] 2 + sin 2 (^) 



(9a) 
(9b) 

(9c) 

• (9d) 



The parameter 77 = ^ is a measure of the Kondo effect 

strength in the lower QD. Eq. (pi) reveals a shift by e Q 
of the center of the resonance away from its bare value 
£u = £f = 0, as a consequence of interference. 



trl(%) 



[e u ~ iT u ] [e d - iT d ] + T u T d cos 2 ( § ) 



• (7) 



In Eq. Q we have set the Fermi level to zero (ep = 0) 
and introduced the coupling strengths T u = 2-^ and 

reduces to the transmission calcu- 
lated by Lopez et alW^ for identical QDs and symmetric 
junctions. 

Having obtained the ingredients to calculate the con- 
ductance in the presence of Coulomb interaction in the 
lower QD, we now consider the regime in which the lower 
arm of the Aharonov-Bohm ring is tuned to the Kondo 
regime and scattering occurs within a continuum (back- 
ground process). In the upper arm transport occurs 
through the discrete QD energy level (resonant process). 
The interference between both processes will give rise to 
a Fano line-shape in the conductance. The transport is 
then characterized by the Fano parameter q, which in 
general is a complex number. The parameter \q\ mea- 
sures the ratio of the resonant scattering to the back- 



Ill. CALCULATIONS & RESULTS 

In this section we analyze the result for the conduc- 
tance [Eq. ^] from the previous section by solving the 
system of equations [Eqs. Q-([6|] numerically. Unless 
stated otherwise we take the wide band limit and the 
symmetric case: T u = T d 
unit of energy. 



jfo, where t Q is taken as the 



A. Kondo channel as continuum scattering path 

We begin by showing that our calculations indicate 
that the Kondo channel through the lower QD is acting 
as the continuum scattering path within the Fano effect. 
Fig. |2|a) shows the different line-shapes for the three 
regimes mentioned in section [TTJ When e d » T d the 
probability of an electron at the Fermi level to tunnel to 
the lower QD is negligible. The conductance is therefore 
determined by electrons transversing the upper QD. This 
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FIG. 2. (a) Normalized conductance, C, for the three different 
regimes: Lorentzian (resonant channel dominates, red dot- 
dashed curve), mixed- valence (resonant channel and Kondo 
channel contribute, purple dashed curve) and Fano-Kondo 
(Kondo channel dominates, black continuous curve), (b)-(e) 
Normalized conductance, G, versus energy levels, e u and ed, 
for different magnetic fluxes through the ring. 



gives a resonance with a Lorentzian shape centered at the 
Fermi level e u = [red point-dashed curve in Fig. |2fa)] 
for an arbitrary magnetic phase shift <fi [Figs. [2jb)-|2fe)]. 
In this regime there is no Kondo effect, rj >> 1, and Eq. 
181) reduces to: 



G 



(¥) 



■n 



(10) 



We also find \q\ » 1, hence the system behaves as if 
there were no lower path in the ring and transport is 
dominated by the resonant channel. 

We now consider the interacting QD to be in the 
mixed- valence regime « — r^. Here, both arms con- 
tribute to the transmission through the structure. The 
interference between the resonant upper path and the 
mixed-valence lower path, acting as a continuum, gives 
rise to a Fano shape [dashed curve in Fig. [2ja)] with 
its symmetry tuned by the magnetic phase shift <j) [Figs. 
[2|b)-[2|e)]. This symmetry is obtain by computing the 
argument of the Fano parameter 9 [Eq. (9d)]: 



arg(g) = arctan 



tan (f) - ^+T cot (f). 



(11) 



In this regime the Fano parameter gives \q\ ~ 1, which 
agrees with the interpretation of \q\ [for the green dashed 
curve in Fig. [2|a) \q\ = 1.2]. 



B. Symmetry change due to Fano-Kondo interplay 

In a transport experiment performed by Verduijn 
et alM using a three-dimensional silicon field effect tran- 
sistor (FinFET), where a few arsenic dopants diffuse to 
the transistor channel, an asymmetric Fano resonance 
for one of the arsenic atoms was observed at zero-bias 
voltage, and speculated that the Fano resonance results 
from the interference between a Kondo transport channel 
and a direct transport process. The geometry consists of 



two parallel paths with an impurity in each one. Us- 
ing parameters similar to the ones in this experiment, 
= 10r w and €d = — 3r^, we obtain qualitative agree- 
ment. Fig. |3|b) shows the calculated conductance for a 
magnetic phase shift of 27r/5, in order to obtain a compa- 
rable degree of symmetry with the reported resonance [in- 
set in Fig. |3|b)]. In this calculation, a value of \q\ = 1.14 
was obtained, which indicates that the measured asym- 
metric conductance profile originates from an underlying 
mixed-valence Kondo path to the upper QD bare reso- 
nance. The degree of symmetry [arg(g)] versus magnetic 
field was calculated in the experiment, and an abrupt 
jump in symmetry at half a period was foundP We rec- 
ognize this sudden symmetry change as a trademark of 
an underlying Kondo channel of the transport process. 
Fig. [3^ a) shows the argument of the calculated Fano 
parameter [Eq. ( fl!] )] for = 10T U . A linear change 
in symmetry occurs when the lower QD is completely 
outside the Kondo regime [black continuous line in Fig. 
|3ja)]. As soon as the QD enters the Kondo regime (even 
in the mixed- valence regime) an abrupt symmetry change 
begins to form at half a period [green dashed line in Fig. 
pta)]. This transition becomes more abrupt as the lower 
QD enters deeper into the Kondo regime (rj << 1)P^ 



(a) 



^0 . 









— €d~- 


=3r d 






---Q= 


=-3r d 













-1 








FIG. 3. Calculations performed for Yd — 10r u . (a) Degree of 
symmetry of the resonance obtained from Eq. (11) for 6d = 
3Fd (black continuous) and €d = —SFd (green dashed), (b) 
Calculated conductance for an energy level of the interacting 
QD as in the experiment (ed = —3Td). Inset in (b) shows the 
reported data and fit of the experiment [J. Verduijn et al., 
Appl. Phys. Lett. 96, 072110 (2010)]. 



In the Fano-Kondo regime the lower QD is tuned to the 
Kondo state by lowering its energy level to « — 
[black continuous curve in Fig. |2ja)]. The shape of the 
conductance then changes from a peak to a dip line- 
shape. As was observed experimentally,^ this conduc- 
tance behavior indicates that the Kondo channel has 
taken over the role of scattering within a continuum in a 
Fano effect, interfering with a discrete scattering process 
mediated by the upper QD. For a magnetic flux away 
from (j) = 27rn, transport is dominated by the Kondo 
channel regardless of the degree of interference [Figs. 
|^b)-|^e)]PESl Within this regime the lower QD is in the 
Kondo state, then A — » — such that 77 << 1. In this 
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limit \q\ = rj « 1 and 



G 



eS + [r u sin 2 (f)]2 



(12) 



In contrast with the constant broadening of the peak 
shape [red point-dashed curve in Fig. [2] (a)] the broad- 
ening of the dip is modulated by sin 2 (f) as seen in Figs. 
[2^b)-[2je). Furthermore, the expression for q [Eq. ( |9d| )] is 
in accordance with the appropriate interpretation of the 
phenomenon [for the black continuous curve in Fig. [2ja) 
kl = 0.3]. 



C. Tuning the electron's path preference 

So far, we have seen that our system can be in three 
different regimes depending on the relative position of the 
energy level of the interacting QD. Our results show that 
the magnetic flux can be used as a tool to tune the con- 
ductance within each of the regimes with characteristic 
landmarks: no effect in the resonant regime, symmetry 
change of Fano resonance in the mixed-valence regime 
and change of the dip broadening in the Fano-Kondo 
regime. It is in this last regime where the magnetic flux 
can have the extra role of tuning the gap between two res- 
onances and influences not only the interference at the 
end of the ring, but affects the probability of the electron 
taking the lower or upper path. 



(a)0 = (b)0 




FIG. 4. Calculations perfomed for ea = — SFd. (a-d) Conduc- 
tance, G, versus the noninteracting QD energy level, e u , for 
different magnetic phase shifts, 0, close to zero, (e) Conduc- 
tance, C, versus e u and magnetic phase 0. Figs, (a)-(d) are 
different horizon tal c uts of (e). (f) Absolute value of the Fano 
parameter [Eq. (|9d[)] for the same parameters as in (e). 



When <p « 2nir the conductance [Eq. ([8|] cannot be 
written as Eq. (12). Even though the conductance for 
small magnetic shifts <\> [Figs. |4^a)-|4jc)] has a distinc- 
tive line-shape, the antiresonance near e u = 0, of width 
r oc e 2 d + f^sin 2 (^), remains present for any value of 
magnetic flux. This magnetic flux independence is due 
to the fact that near e u = the upper discrete path is 
close to resonance, such that the transport characteris- 
tics through it remain robust under the influence of a 
magnetic flux. Therefore, the only effect remaining is 



the interference with the Kondo channel, giving the al- 
ready mentioned dip line-shape [Figs. [4ja)-|4jc)]. More- 
over, two sharp resonances appear when <j) w 2nir [Figs. 
|4^a)-j4^c)] . Since these two resonances originate from the 
contribution of the discrete energy level path and the 
Kondo path [\q\ « 1 in Fig. |4jf) as indicated by the con- 
tour dashed line mapping both resonances] , we speculate 
that they appear because of a hybridization of the non- 
interacting QD state with the Kondo state, which leads 
to a gap opening between two hybrid states. This gapped 
system could lead to interesting physics as we control the 
state of the system (under or above the gap) as well as 
the size of the gap, by tuning the non- interacting QD 
energy level and the magnetic flux, respectively. 

We can also use the magnetic field to alter the elec- 
tron's path preference, even if the interacting QD is deep 
in the Kondo regime and this path becomes the pre- 
ferred channel for the scattering process. Fig. |4je) shows 
how the conductance behavior changes near </> = and 
the corresponding absolute value of the Fano parame- 
ter [Fig. [4^f)] indicates that the transport mechanism 
changes from dominance of the Kondo channel (\q\ « 0) 
to a mechanism in which both scattering paths contribute 
(\q\ « 1). The effect is a consequence of the decrease 
of the Kondo temperature via the manipulation of the 
Kondo- QD tunnel coupling by the magnetic flux. The 
tunnel coupling decreases until it becomes equally prob- 
able to tunnel to either QD, such that both scattering 
paths contribute significantly to the conductance. The 
dependence of the Kondo temperature on magnetic flux 
can be derived by writing the Fano form for conductance 
[Eq. ^] in terms of the interacting energy level e^, such 
that the expression found will describe the Kondo reso- 
nance. Then, the derived state broadening (r) is sub- 
stituted into the expression for the Kondo temperature 
using SBMF HamiltonianpSl 



B-l-K 



~2 



f3 



sin 



C 2 + i C 2 + i 



(13) 



with C — fr"- , k>B the Boltzman constant and Tk the 
Kondo temperature. 



IV. CONCLUSIONS 

In this paper, we have studied coherent transport 
through an Aharonov-Bohm ring structure with two em- 
bedded QDs. In particular, we have analyzed the in- 
fluence of a Coulomb interacting channel to a discrete 
scattering transport process in an Aharonov-Bohm con- 
figuration. Three regimes were identified, depending on 
the contribution of the Kondo channel to the conduc- 
tance. These regimes yield different conductance line- 
shapes: Lorentzian, asymmetric Fano and symmetric dip. 
The influence of a magnetic phase shift to these regimes 
is to tune the Fano resonance asymmetry and to tune 
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the broadening of the symmetrical dip. We have shown 
that the Fano resonance line-shape can be explained by 
an underlying scattering mechanism through a Kondo 
channel, and that the shape can vary from an asymmet- 
ric resonance to a symmetric dip as we drive the Coulomb 
blocked QD further into the Kondo regime. We have also 
shown that we can identify the presence of an underlying 
Kondo channel by an abrupt change in symmetry of the 
discrete channel resonance. Furthermore, we have found 
agreement with experimental data that support our find- 
ings. Finally, we predict a gap opening between two res- 
onances that we speculate is a consequence of gapped 
hybrid states. We also predict that, even if the interact- 
ing channel is fully in the Kondo regime, we can use the 
magnetic flux to diminish its contribution by lowering the 
characteristic Kondo temperature (Kondo state broaden- 



ing), producing an alteration in the electron's path pref- 
erence. 
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